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The ground-state and low-energy excitations of quantum Hall systems are studied by the 
density matrix renormalization group (DMRG) method. From the ground-state pair corre- 
lation functions and low-energy excitions, the ground-state phase diagram is determined, 
which consists of incompressible liquid states, Fermi liquid type compressible liquid states, 
and many kinds of CDW states called stripe, bubble and Wigner crystal. The spin transition 
and the domain formation are studied at z/ = 2/3. The evolution from composite fermion 
liquid state to an excitonic state in bilayer systems is investigated at total filling factor u = 1. 

§1. Introduction 

In two dimensional systems, applying perpendicular magnetic field strongly mod- 
ifies the wave function of electrons leading to many interesting phenomena at low 
temperatures. The fractional quantum Hall effects^' are typical example, where in- 
compressible ground states are realized only at some fractional fillings of Landau 
levels.'^''^ Since fractional quantum Hall effects are observed only in high quality 
samples, the Coulomb interaction between the electrons is thought to be essential 
rather than random potentials from impurities. This is contrasted with the case of 
integer quantum Hall effect where random potentials are essential. The importance 
of the Coulomb interaction in high magnetic field is followed from the increase in 
the energy scale of the Coulomb interaction. The wave function is scaled by the 
magnetic length £ = y^h/eB, which is equivalent to the classical cyclotron radius Tc 
in the lowest Landau level. The increase in the magnetic field decreases the mag- 
netic length and enhances the energy scale of the Coulomb interaction between the 
electrons, e^/ei. 

At typical magnetic field of lOT, £ is about 8nm, which is still much larger 
than the atomic length of O.lnm. Since the conduction electrons are on the positive 
background charge from ions over length scale of £, the positive charge may be 
simplified to be uniform. Then the system is equivalent to the electron gas in a 
magnetic field, and £ becomes unique length scale of the system. 

In the magnetic field, the kinetic energy is scaled by the cyclotron frequency 
ujc = eB/ni, which is determined by the magnetic field B. The quantization of the 
wave function discretizes the classical cyclotron radius Tc, which also discretizes the 
kinetic energy and makes Landau levels E = hiOc{n + 1/2). This means that the 
macroscopic number of electrons have the same energy in each Landau level, and 
large-scale degeneracy appears in the ground state. This macroscopic degeneracy is 
lifted by the Coulomb interaction between the electrons and various types of liquid 
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stated '13 '0 'El and CDW stateP'^i^'IS'lSJ are realized depending on the filling of 
the Landau levels. 

Since the ground state has macroscopic degeneracy in the limit of weak Coulomb 
interaction, standard perturbation theories are not useful. Thus numerical diagonal- 
izations of the many body Hamiltonian have been used to study this system. Since 
numerical representation of the Hamiltonian needs complete set of many body basis 
states, we divide the system into unit cells with finite number of electrons in each 
cell. The properties of the infinite system are obtained by the finite size scalings. 
However, the number of many body basis states increases exponentially with the 
number of electrons. For example, when we study the ground state at v = 1/3 
with 18 electrons, each unit cell has 54 degenerated orbitals. The number of many 
body basis states is given by the combination of occupied and unoccupied orbitals, 
54C18 ~ 10^^, which is practically impossible to manage by using standard numerical 
method such as exact diagonalization. 

To study systems with typically more than 10 electrons, we need to reduce the 
number of many body basis states. For this purpose, we use the density matrix renor- 
malization group (DMRG) method, which was originally developed by S. White in 
1992.II3:[il This method is a kind of variational method combined with a real space 
renormalization group method, which enables us to obtain the ground-state wave 
function of large-size systems with controlled high accuracy within a restricted num- 
ber of many body basis states. The DMRG method has excellent features compared 
with other standard numerical methods. In contrast to the quantum Monte Carlo 
method, the DMRG method is free from statistical errors and the negative sign prob- 
lem, which inhibit convergence of physical quantities at low temperatures. Compared 
with the exact diagonalization method, the DMRG method has no limitation in the 
size of system. The error in the DMRG calculation comes from restrictions of the 
number of basis states, which is systematically controlled by the density matrix cal- 
culated from the ground-state wave function, and the obtained results are easily 
improved by increasing the number of basis states retained in the system. 

The application of the DMRG method to two-dimensional quantum systems 
is a challenging subject and many algorithms have been proposed. Most of them 
use mappings on to effective one-dimensional models with long-range interactions. 
However, the mapping from two-dimensional systems to one-dimensional effective 
models is not unique and proper mapping is necessary to keep high accuracy. In two- 
dimensional systems under a perpendicular magnetic field, all the one-particle wave 
functions ^nx{x, y) are identified by the Landau level index and the x-component 
of the guiding center, X, in Landau gage. The guiding center is essentially the center 
coordinate of the cyclotron motion of the electron and it is natural to use X as a 
one-dimensional index of the effective model. More importantly, X is discretized in 
finite unit cell of x Ly through the relation to y-momentum, X = kyl^ , which 
is discretized under the periodic boundary condition, ky = 2Tm/Ly with n being 
an integer. Therefore, the two-dimensional continuous systems in magnetic field are 
naturally mapped on to effective one-dimensional lattice models, and we can apply 
the standard DMRG method.^5,1 

This method was first applied to interacting electron systems in a high Landau 
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level and the ground-state phase diagram, which consists of various CDW states 
called stripe, bubble and Wigner crystal, has been determined !^''^ The ground 
state and low energy excitations in the lowest and the second lowest Landau levels 
have also been studied by the DMRG and the existence of various quantum liquid 
states such as Laughlin state and charge ordered states called Wigner crystal have 
been confirmed and new stripe state has been proposed.'^''^ 

In the following, we first explain the effective one-dimensional Hamiltonian used 
in the above studies and then show the results obtained for the spin polarized single 
layer system. We next review recent study on the spin transition and domain forma- 
tion at u = 2/3,'22J and finally explain the results on bilayer quantum Hall systems 
at u = Ij'^D where crossover from a Fermi liquid state to an excitonic incompressible 
state occurs. 



§2. DMRG method 



Here we briefly describe how the effective ID Hamiltonian is obtained from 2D 
quantum Hall systems.!^ To describe the many body Hamiltonian for a interacting 
system, we first need to define one-particle basis states. Here, we use the eigenstates 
of free electrons in a magnetic field as one-particle basis states and represent the 
wave function 'Fiyx{x,y) in Landau gauge: 



^Nx{x,y) = Catcxp 



ikyy 



{x - Xf 
2£^ 





'x - X' 




I 



(2-1) 



where are Hermite polynomials and Cm is the normalization constant. Then 
all the eigenstates 'l^Nx{x,y) are specified using two independent parameters and 
X; N is the Landau level index and X is the x-component of the guiding center 
coordinates of the electron. Since the guiding center X is related to the momentum 



ky as X^ 



and ky is discretized under the periodic boundary conditions, the 



guiding center X takes only discrete values 

Xn = 2TTfn/Ly, 



(2-2) 



where Ly is the length of the unit cell in the y-direction. 

If we fix the Landau level index A^, all the one-particle states are specified by one- 
dimensional discrete parameter X^- Since many body basis states are product states 
of one-particle states, they are also described by the combinations of Xn of electrons 
in the system. Thus the system can be mapped on to an effective one-dimensional 
lattice model. 

The macroscopic degeneracy in the ground state of free electrons in partially 
filled Landau level is lifted by the Coulomb interaction 



V{r) 



e 
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(2-3) 



The Coulomb interaction makes correlations between the electrons and stabilizes 
various types of ground states depending on the filling v of Landau levels. When 
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Fig. 1. Schematic diagrams for (a) infinite system algoritfim and (b) finite system algoritlim of the 
DMRG method. • represents a one-particle orbital in a given Landau level. Bl and B_r are left 
and right blocks, respectively. 



the magnetic field is strong enough so that the Landau level splitting is sufficiently 
large compared with the typical Coulomb interaction e^/(e£), the electrons in fully 
occupied Landau levels are inert and the ground state is determined only by the 
electrons in the top most partially filled Landau level. 
The Hamiltonian is then written by 

H = S cJjCn + — ^nin2n3n4,CniCl^2CngCn4^, (2'4) 

n n\ n2 ris rn 

where we have imposed periodic boundary conditions in both x- and y-directions, 
and S is the classical Coulomb energy of Wigner crystal with a rectangular unit cell 
of Lx X Ly ^ cli is the creation operator of the electron represented by the wave 
function defined in equation (j2-ip with X = Xn- ^nin2n3n4 are the matrix elements 
of the Coulomb interaction defined by 

A -A' ^""^^ 



J J jL^ ni-ni,qyLy/2'iT 



[L^{q^f/2)]\^p 



^[n^ — n-i) 

2 ^ M 



(2-5) 



where Ln{x) are Laguerre polynomials with being the Landau level index.'^ 
^ni,n2 ~ 1 when ni = n2(mod M) with M being the number of one-particle states 
in the unit cell, which is given by the area of the unit cell 2ttM£^ = L^Ly. 

In order to obtain the ground-state wave function we apply the DMRG method.!^ 
As shown in Fig. [1] (a), we start from a small-size system consisting of only four 
one-particle orbitals whose indices n are 1, 2, M — 1, and M, and we calculate the 
ground-state wave function. We then construct the left block containing one-particle 
orbitals of n = 1 and 2, and the right block containing n = M — 1 and M by using 
eigenvectors of the density matrices which are calculated from the ground-state wave 
function. We then add two one-particle orbitals n = 3 and M — 2 between the two 
blocks and repeat the above procedure until M one-particle orbitals are included 
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Fig. 2. Eigenvalues Wa of the density matrix for two-dimensional system of 54 orbitals with 18 
electrons. Sum of Wa is equivalent to the norm of the ground-state wave function and normalized 
to be unity. 



in the system. We then apply the finite system algorithm of the DMRG shown in 
Fig- [I] (b) to refine the ground-state wave function. After we have obtained the 
convergence, we calculate correlation functions to identify the ground state. 

The ground-state pair correlation function g{r) in guiding center coordinates is 
defined by 



9{r) 



LxLy 



1 



.{^\Y^d{r-R, + R,)\'F) 



(2-6) 



where Ri is the guiding center coordinate of the ith electron, and it is calculated 
from the following equation 
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(2-7) 



where ^ is the ground state and A'^e is the total number of electrons. 

The accuracy of the results depends on the distribution of eigenvalues of the 
density matrix. A typical example of the eigenvalues of the density matrix for system 
of M = 54 with 18 electrons is shown in Fig. [21 which shows an exponential decrease 
of eigenvalues Wa- In this case accuracy of 10~^ is obtained by keeping more than 
one hundred states in each block. 

§3. Single layer system 



Here we present diverse ground states obtained by the DMRG method applied 
to the single layer quantum Hall systems. In the limit of strong magnetic field, 
the electrons occupy only the lowest Landau level N = 0. In this limit, fractional 
quantum Hall effect (FQHE) has been observed at various fractional fillings.'^ The 
FQHE state is characterized by incompressible liquid with a finite excitation gap.'^ 
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Fig. 3. The lowest excitation gap at various v in the lowest Landau level. Relatively large excitation 
gap is obtained at fractional fillings u = n/{2n + 1). The excitation gap is in units of e'^/{e£). 

These FQHE states are confirmed by the DMRG calculations, where relatively 
large excitation gaps are obtained at various fillings between i/ = 1/2 and S/lcP 
as shown in Fig. [3l We clearly find large excitation gaps at fractional fillings v = 
1/3,2/5,3/7,4/9 and 5/11, which correspond to primary series of the FQHE at 
I' = n/{2n + 1). The pair correlation function at = 1/3 is presented in Fig. U 
which shows a circularly symmetric correlation consistent with the Laughlin's wave 
function.^ 

In the limit of low filling ^ 0, mean separation between the electrons becomes 
much longer than the typical length-scale of the one-particle wave function. In this 
limit the quantum fluctuations are not important and electrons behave as classical 
point charges. The ground state is then expected to be the Wigner crystal. The 
formation of the Wigner crystal is also confirmed by the DMRG calculations at low 
fillings as shown in Fig. [5] (a). The i^-dependence of the low energy spectrum shows 
that the first-order transition to Wigner crystal occurs at z> ~ 1/7.^^ ' 

With decreasing magnetic field, electrons occupy higher Landau levels. In high 
Landau levels, the one-particle wave function extends over space leading to effective 
long range exchange interactions between the electrons. The long range interaction 
stabilizes CDW ground states and various types of CDW states called stripe and 
bubble are predicted by Hartree-Fock theory."^ These CDW states are confirmed by 
the DMRG calculations as shown in Figs. [5] (b) and (c), where two-electron bubble 
state and stripe state are obtained at = 8/27 and 3/7, respectively, in the N = 2 
Landau level. Although the CDW structures are similar to those obtained in the 
Hartree-Fock calculations, the ground state energy and the phase diagram are sig- 
nificantly different. The DMRG results are consistent with recent experiments,'^^ 
and the discrepancy is due to the quantum fluctuations neglected in the Hartree-Fock 
calculations. 

The ground-state phase diagram obtained by the DMRG is shown in Fig. [6l In 
the lowest Landau level, we find many liquid states at fractional fillings and around 
z^ = 1/2. Nevertheless, CDW states dominate over the whole range of filling in 
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higher Landau levels. This difference in the ground state phase diagram comes from 
different effective interactions between the electrons. In the lowest Landau level, 
the one particle wave function is localized within the magnetic length ^, that yields 
strong short-range repulsion between the electrons. Since quantum liquid states such 
as Laughlin state are stabilized by the strong short-range repulsion, liquid states are 
realized in the lowest Landau level. In higher Landau levels, however, the wave 
function extends over space with the increase in the classical cyclotron radius re- 
Thus the short-range repulsion is reduced and liquid states become unstable. As 
shown in Fig. [7] (a), the real space effective interaction between the electrons in 
higher Landau levels has a shoulder structure around the distance twice the classi- 
cal cyclotron radius. This structure of effective interaction makes minimum in the 
Coulomb potential near the guiding center of the electron as shown in Fig. [7| (b) and 
stabilizes the clustering of electrons. This is the reason why stripe and bubble states 
are realized in higher Landau levels .^^^ 



Laughlin state 
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Fig. 4. Pair correlation function g{r) a.t p = m the lowest Landau level. The length is in units 
oil. 



(b) two-electron 
(a) Wigner crystal bubble 




Fig. 5. Pair correlation functions g{v) in guiding center coordinates, (a) Wigner crystal realized in 
an excited state at = 1/6 in the lowest Landau level. The number of electrons in the unit 
cell iVe is 12. (b) Two-electron bubble state a.t v = 8/27 in = 2 Landau level. = 16. (c) 
Stripe state at u = ?>/7 in N = 2 Landau level. Ne = 18. 



8 



Naokazu Shibata 



cr.oc 2/7 FQHE FQHE,,^' 

AT—Ct FQHE . ■ TT 4/9 

^-tJ , , Stnpell 1/3 2/5 3/7 

Wigner crystal ''J 



_L 



liquid 





N=l 



0.1 



0.4 



0.5 



0.2 0.3 
V 

FQHE FQHE 
Wigner crystal !« stripe H 1/3 stripe I paWng 

' -h — u 

0.1 0.2 0.3 0.4 0.5 



_l l_ 

0.2 0.3 
V/v 



Wigner crystal 



2-electron 
bubble 



stripe I 



0.1 



0.2 0.3 



0.4 



0.5 



Fig. 6. The ground state phase diagram obtained by the DMRG method. A'^ is the Landau level 
index and in the filling factor of the A''th Landau level. 
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Fig. 7. (a) Effective interaction between the electrons in the A'^th Landau level. Rc is the classical 
cyclotron radius, (b) Coulomb potential made by two electrons separated by Ax. 



§4. Spin transitions 



In two dimensional systems, strong perpendicular magnetic field completely 
quenches the kinetic energy of electrons. Since the kinetic energy is independent 
of the spin polarization, the exchange Coulomb interaction easily aligns the electron 
spin. The ferromagnetic ground state at u = 1/q {q odd) is thus realized even in 
the absence of the Zeeman splitting.!^ At the filling u = 2/3 and 2/5, however, the 
paramagnetic ground states compete with the ferromagnetic state, and the Zeeman 
splitting Az = gfJ-sB induces a spin transition.!^ Such a spin transition in fractional 
quantum Hall states has been naively explained by the composite fermion theory.l^S} 
The composite fermions are electrons coupled with even number of fluxes. These 
fluxes effectively reduces external magnetic field and the = p/{2p it 1) fractional 
quantum Hall effect (FQHE) state is mapped on to the v' = p integer QHE state 
of composite fermions. The spin transitions at u = 2/3 and 2/d^ correspond to 
the spin transition at = 2, where the Zeeman splitting corresponds to the effective 
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Fig. 8. Lowest energies for fixed polarization ratio P as a function of magnetic field B at filling 
factor V — 2/3 in units of /{et). The total number of electron is 20. The aspect ratio is fixed 
at 2.0. The g-factor is 0.44. 
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Fig. 9. Charge gap of = 2/3 spin polarized states (□), unpolarized states (•), and partially 
polarized states (o) for various Ne and aspect ratios L^jLy. Ac is in units of /{et). 



Landau level separation, and the energy levels of the minority spin state in the lowest 
Landau level and the majority spin state in the second lowest Landau level coincide. 
Extensive experimental2it'|23.[23,l3l,ED,[3l,|33,[3l,[35j theoreticaPS'EBS® 

studies have been made on this transition. Nevertheless, there is no clear theoretical 
consensus on this issue. This is due to the difficulties of studies in this system. A 
number of states possibly compete in energy, and large enough systems are needed to 
see non-uniform structures in the partially polarized states. Here we use the DMRG 
method,'!^!' and study the spin transition and the spin structures in large system to 
clarify the nature of the spin transition at = 2/3. 

We first calculate the energy at various polarization P as a function of the 
Zeeman splitting, Az = gfJ-B. The obtained results are shown in Fig. [8l In the 
absence of the Zeeman splitting, the unpolarized state (P = 0) is the lowest. The 
energy of the polarized state (P > 0) monotonically increases as P increases. With 
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the increase in Zeeman splitting A^, however, the energy of polarized state decreases 
and the fully polarized state {P = 1) becomes the lowest. Figure [8] shows that the 
transition from the unpolarized state to the fully polarized state occurs at ~ 6T 
which is roughly consistent to the earlier work done in a spherical geometry.^ In 
the present calculation on a torus, all partially polarized states (0 < P < 1) are 
higher in energy than the ground states (P = or 1). This feature is independent of 
the size of the system and the aspect ratio L^/Ly, and indicates phase separations 
of P = and P = 1 in partially polarized states. 

The unpolarized state of P = and the fully polarized state of P = 1 are both 
quantum Hall states with finite charge excitation gap, which is defined by 

Z\e(P) = E{N^ + 1,P) + E{N^ - 1,P) - 2E{N^,P), (4-1) 

where is the number of one-particle states in the lowest Landau level. The filling 
factor u is then given by N^/N^. The charge gap for various N^j, and aspect 
ratios of the unit cell is presented in Fig. [9j In this figure, the gap Ac seems to 
vanish for partially polarized state P ~ 1/2 in the limit of N^^ — > oo. This result 
clearly indicates that partially polarized state with P ~ 1/2 is a compressible state 
in contrast to the incompressible states at P = and 1, where Ac remains to be 
finite in the limit oi Nc ^ oo. 

To study the spin structure in the partially polarized states, we next calculate 
the pair-correlation function defined by 

where a = ±1/2 is the spin index and No- is the number of electrons with spin 
a. The spin structures in partially spin polarized states are clearly shown in the 
pair correlation function between minority spins. Namely, if unpolarized regions 
are formed in the partially polarized states, then electrons with minority spins are 
concentrated in the unpolarized regions. This concentration of the minority spins is 
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shown in Fig. IIOI which shows gii{x^y) for partiahy polarized states at (a) P = 0.8, 
(b) P = 0.6, (c) P = 0.5, and (d) P = 0.4. When P is close to 1, for example P = 0.8 
shown in Fig. [TOlfa). a pair of minority spins is found only near the origin. As the 
polarization ratio P decreases, minority spins make a domain around the origin, and 
two domain walls along the y-direction are formed. These domain walls move along 
x-direction and the domain of minority spin finally covers entire unit cell in the limit 
of P = 0. This change in the size of the domain is consistent with the expectation 
that the domain in Fig. [10] corresponds to the unpolarized spin singlet region where 
the density of up-spin electrons and the down-spin electrons are the same. 

To confirm the separation of the unpolarized and polarized spin regions, we 
next consider the local electron density of up-spin electrons i''\{x) and down-spin 
electrons i^iix). Figure [TT] shows i']{x) and ^[{x) for partially polarized states with 
P = 0.2, 0.4, 0.6 and 0.7. Here i"\{x) and vi{x) are scaled to be the local filling 
factor of the lowest LL. Thus, the total local electron density v^{x) + viix) is almost 
2/3. In this figure the separation to two regions is clearly seen; the unpolarized spin 
region around where both v-^ and V]^ are close to 1/3, and the fully polarized 

spin region around x ~ or equivalently x ~ L^,, where is almost 2/3 while 
is close to 0. These results confirm the separation of the unpolarized and polarized 
spin regions as expected from the pair correlation functions shown in Fig. [TTl 

The polarized and unpolarized spin regions are separated by the domain walls 
whose width is about M. This means that the phase separation is realized only 
for systems whose size of the unit cell Lx, {Ly) is larger than twice the width of 
domain wall; Lx,{Ly) > 8i. Indeed, exact diagonalization studies up to A^'e = 8 
electrons have never found the phase separation at v = 2/3.® We have found the 
phase separation only for large systems with Nf^ > 12. We note that above behavior 
is generic over the aspect ratio. In an ideal system, the two states separate into 
two regions even when the system size is infinitely large. In experimental situations, 
however, multi-domain structures are realized due to the inhomogeneity and coupling 
with randomly distributed nuclear spins. 
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The DMRG study on the ground state energy for various polarization P shows 
that the ground state at = 2/3 evolves discontinuously from the unpolarized P = 
state to the fully polarized P = 1 state as the Zeeman splitting increases. In partially 
polarized states < P < 1, the electronic system separates spontaneously into two 
states; the P = and the P = 1 states. These two states are separated by the 
domain wall of width 4£. Since the energy of the domain wall is positive, the partially 
polarized states always have higher energy than that of -P = 1 or P = states. We 
think this is the reason of the direct first order transition from P = to P = 1 state 
in the ground state. 

It is useful to compare our result with the spin transition at u = 2 which occurs 
when minority spin states in the lowest LL and majority spin states in the second 
lowest LL cross by varying the ratio of the Zeeman and Coulomb energy. The 
ground state at = 2 is thus a fully polarized state or a spin singlet state. In 
analogous to the case oi u = 2/3 the transition between them is first order 
and spin domain walls have been found in high energy states.'^D' This analogy can 
be expected, because the v = 2 states and the v = 2/3 states are connected in the 
composite fermion theory,'^''^ although the effective interaction between composite 
fermions is different from that for electrons. 

§5. Bilayer system 

The properties of quantum Hall systems sensitively depend on the magnetic field, 
and various types of ground states including incompressible liquids,'^ compressible 
liquids,H3'S3 spin singlet liquid, CDW states called stripes, bubbles, and Wigner 
crystal are realized depending on the filling of Landau levels. In bilayer quantum 
Hall systems, additional length scale of the layer distance d, and the degrees of 
freedom of layers make the ground state much more diverse and interesting.^ 

Excitonic phase, namely Haplerin's l^'i^i.i state, is one of the ground states real- 
ized in bilayer quantum Hall systems at total filling u = 1 at small layer separation 
d, where electrons and holes in different layers are bound with each other due to 
strong interlayer Coulomb interaction. This excitonic state has recently attracted 
much attention because a dramatic enhancement of zero bias tunneling conductance 
between the two layers,H9 and the vanishing of the Hall counterflow resistance are 
observed .SSJ '33 As the layer separation d is increased, the excitonic phase vanishes, 
and at large enough separation, composite-fermion Fermi-liquid state is realized in 
each layer. 

Several scenarios have been proposed for the transition of the ground state as 
the layer separation incr eases. H51''H21''l5S,[5lJ,|52),[53li,|54li jjowever how the excitonic state 

develops into independent Fermi-liquid state has not been fully understood. In this 
section we investigate the ground state of z/ = 1 bilayer quantum Hall systems 
by using the DMRG method.^ We calculate energy gap, two-particle correlation 
function g(r) and excitonic correlation function for various values of layer separation 
d, and show the evolution of the ground state with increasing d. 
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Fig. 12. The exciton correlation ol bilayer quantum Hall systems &t v = 1. The solid line represents 
3gx(M/2). The dashed line represents gex(M/2 — 1). 



The Hamiltonian of the bilayer quantum Hah systems is written by 

i<j q 
i<j q 



(5-1) 



where Ri^i are the two-dimensional guiding center coordinates of the ith electron 
in the layer-1 and -R2,j are that in the layer-2. The guiding center coordinates 
satisfy the commutation relation, [Rj,R^] = li'^Sjk- V{q) = 27re^/(e(7) is the Fourier 
transform of the Coulomb interaction and the wave function is projected on to the 
lowest Landau level. We consider uniform positive background charge to cancel 
the component at q = 0. We will assume zero interlayer tunneling and fully spin 
polarized ground state. 

In the limit of d = 0, electrons in different layers can not occupy the same 
position because of the strong interlayer Coulomb repulsion. The strong interlayer 
repulsion makes electron-hole pairs, which is called excitons whose degrees of freedom 
are represented by interlayer dipoles or pseudo-spins at total filling = 1. The 
Coulomb exchange interaction aligns the interlayer dipoles (the pseudo-spins) leading 
to the macroscopic coherence of the excitons and Haplerin's state is realized. 

To confirm the coherence of the excitons, we calculate the exciton correlation 
defined by 

5ex(n) = ^j^^/ (^|ci«C2,«4.oCi.o|j^), (5-2) 



where \^) is the ground state and c\ „ „) is the creation operator of the electrons 
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Fig. 13. The pseudo-spin excitation gap Zips of bilayer quantum Hall systems at the total filling 
factor u — 1. The dashed lines are guide for the eye. 

in the nth one-particle state defined by 



in the layer-1 (layer-2). Xn = nLx/M = kyi'^ is the x-component of the guiding 
center coordinates and is the length of the unit cell in the x direction. M is the 
number of one-particle state in each layer. A^i and N2 are the number of electrons 
in the layer-1 and layer-2, respectively, and we impose the symmetric condition 
of A'^i = A''2. Since gcx{n) represents the correlations between the two excitons at 
X = and X = Xn, lim„^oo dcxin) 7^ indicates existence of macroscopic coherence 
of excitons. 

As is shown in Fig. [T2l gc^{n) tends to 1 as d ^ 0, that confirms the macroscopic 
coherence of excitons at d = 0. Indeed, Haplerin's state has the macroscopic 

coherence of excitons and gcx{n-) = 1 independent of n. In this figure we have 
shown gexiM/2) instead of lim„^oo 5ex (^t-) > because the largest distance between the 
two excitons is Lx/2 in the finite unit cell of x Ly under the periodic boundary 
conditions. In order to check the size effect, we also plot gcxiM/2 — 1) with the 
dashed line. Since the difference between gc^{M/2) and gcxiM/2 — 1) is small, we 
expect gex{M/2) well represents the macroscopic coherence in the limit of — > 00. 
With increasing d/£, the excitonic correlation decreases monotonically and finally 
falls down to negligible value at d/i ~ 1.6. 

The presence of macroscopic coherence of excitons shown in the Fig.ll2lmeans the 
existence of ferromagnetic order of the interlayer dipoles (the pseudo-spins). Since 
interaction between the interlayer dipoles has SU(2) symmetry at d = 0, corrective 
gapless excitations called pseudo-spin waves are expected. Even in the case of finite 
layer distance d, the Hamiltonian has continuous XY symmetry, and gapless pseudo- 
spin wave excitations are still expected. This is confirmed by the size dependence of 
the pseudo-spin excitation gap shown in Fig. [131 where the pseudo-spin excitation 




(5-3) 
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Fig. 14. The charge excitation gap Ac of bilayer quantum Hall systems at the total filling factor 
V = 1. The dashed lines are guide for the eye. 



gap Aps = E{Ni + 1,N2 — — E{Ni, N2, M) in finite system decreases as a 

function of l/L^- 

In contrast to the pseudo-spin excitation gap Aps, the charge excitation gap Ac 

defined by A^ = E{Ni,N2,M - 1) + E{Ni,N2, M + 1) - 2E{Ni,N2, M) seems to 

be finite even in the limit of L^; ^ cxd for small d as shown in Fig. [TH The charge 

(1 2) 

excitation brakes at least one electron-hole pair and it needs energy of order Vq ' 
which is the pseudopotential between the electrons in different layers whose relative 
angular momentum is 0. This pseudopotential decreases with the increase in the 
layer distance d, and thus the charge gap decreases with the increase in d. 

We next see the lowest excitation gap in a fixed size of system. Figure [15] shows 
the result for Nf> = 24. The aspect ratio L^/Ly = 1.8 is chosen from the minimum 
of the ground state energy with respect to Lx/Ly around d/l = 1.8, where minimum 
structure appears in the ground state energy. 

We can see clear excitation gap of finite system for d/H. < 1.2, where excitonic 
ground state is expected both theoretically and experimentally.lSJ.I13E7),48),49),50), 53,lll,IS3,[Sl) 
The excitation gap rapidly decreases with increasing d/i from 1.2, and it becomes 
very small for d/£ > 1.6. This behavior is consistent with experiments.'^ Although 
the excitation gap for d/i > 1.7 is not presented in the figure for A^e = 24 because 
of the difficulty of the calculation of excited states in large system, we do not find 
any sign of level crossing in the ground state up to d/i ~ 4, where two layers are 
almost independent. These results suggest that the excitonic state at small d/i con- 
tinuously crossovers to compressible state at large d/i, that is consistent with the 
behaviors of exciton correlations gc^{M/2) in Fig. 1121 which shows gc^{M/2) contin- 
uously approaches zero around d/i ~ 1.6. In the present calculation it is difficult to 
conclude whether the gap closes at finite d/i ~ 1.6 in the thermodynamic limit or 
excitonic state survives with exponentially small finite gap even for large d/i > 1.6. 
We have calculated the excitation gap in different size of systems and aspect ratios, 
and obtained similar results as shown in the inset of Fig. [T5l 
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Concerning the first excited state, however, Fig. [T5l shows a level crossing at 
d/£ ~ 1.2, where we can see sudden decrease in the excitation gap. We expect that 
the lowest excitation at d/i < 1.2 is the pseudo-spin excitation whose energy gap 
decreases with the increase in the size of system and tends to zero in the limit of 
large system. On the other hand, the lowest excitation at d/i > 1.2 shown in Fig. [15] 
is expected to be the excitation to the roton minimum which corresponds to the 
bound state of quasiparticle and quasihole excitatins, whose energy increases with 
the decrease in d/i. This change in the character of the low energy excitations at 
d/£ ~ 1.2 will be confirmed in a clear change in correlation functions in the excited 
state as shown later. We note that the position of the level crossing in the first 
excited state itself depends on the size of system because the pseudo-spin excitation 
gap decreases with the increase in the system size. However, the change in the 
character of the low energy excitations of finite systems is expected to remain even 
in the limit of large system, since the spectrum weight of pseudo-spin waves transfers 
to high energy with the increase in d. 

We next calculate pair correlation functions of the electrons to see detailed evo- 
lution of the ground-state wave function. The interlayer pair correlation functions 
are defined by 

<?i2(r) ^ ^{^\ E '^(^ + ^1.- - R2,m)m^ (5-4) 

^ n m 

where is the ground state. We present Agi2{r) in Fig. 1161 which is defined by 

^5i2(r) = I {guir') - m\r'\ - r) dr' , (5-5) 

where r' is the two-dimensional position vector in each layer. Agi2{r) represents the 
difference from the uniform correlation of independent electrons. 
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Fig. 16. The inter-layer pair correlation function of electrons in the ground state of bilayer quantum 
Hall systems at v = 1. Ne — 24 and Ly/Lx ~ 1.6. 



At d/i = we find clear negative Agi2i'>') around r/£ = 1, which is a characteris- 
tic feature of the excitonic state made by the binding of electrons and holes between 
the two layers. The binding of one hole means the exclusion of one electron caused 
by the strong interlayer Coulomb repulsion. The increase in the layer separation 
weakens Coulomb repulsion between the two layers and reduces |Z\5i2(r)| around 
r/£ = 1. 

The decrease in the interlayer correlation |Z\5(i2(r)| opens space to enlarge corre- 
lation hole in the same layer and reduce the Coulomb energy between the electrons 
within the layer. This is shown in Fig. 1171 which shows the pair correlation functions 
of the electrons in the same layer defined by 

511 W = ^ C^l Yl + ^1.- - Ri,mW)^ (5-6) 



^9n{r) = J (511 (r-O - m\r'\ - r) dr' . (5-7) 

The obtained results indeed show that the correlation hole in the same layer around 
r/i ^ 1 is enhanced with the increase in d/i contrary to the decrease in size of 
interlayer correlation hole in Fig. [TBI The correlation hole in the same layer mono- 
tonically increases in size up to d/i = 1.8, and then it becomes almost constant. 
The correlation function gii{r) for d/£ > 1.8 is almost the same to that of = 1/2 
monolayer quantum Hall systems realized in the limit of d/i = oo. This is consistent 
with the almost vanishing excitation gap and exciton correlation at d/i > 1.8 shown 
in Figs. [Eland [13 

Figure [T7I also shows that the growing of the correlation hole around r/i ~ 1.5 
is accompanied with the increase in Agii{r) around r/i ~ 4. The distance Ai is 
comparable to the approximate mean distance between the electrons 3.54£ estimated 
from {L^Ly/Nif/^ = {lirL/NiY/^i. This means the electrons in the same layer 
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Fig. 17. The intra-layer pair correlation lunction ol electrons in the ground state of bilayer quantum 
Hall systems &t v = 1. Ne — 24 and Ly/Lx — 1.6. 

tend to keep distance of about 4^ from other electrons with the large correlation hole 
around r/£ ~ 1.5 for d/i ^ 1. This is consistent with the formation of composite 
fermions at d = oo, where two magnetic flux quanta are attached to each electron, 
which is equivalent to enhance the correlation hole around each electron in the same 
layer to keep distance from other electrons. 

The large correlation hole in gii{r) attracts electrons in the other layer as shown 
in Fig. [TBI where we find a clear peak in Agi2{r) at r/£ 3. This peak at r/i ~ 3 
is comparable to the neighboring correlation hole at r/i ~ 1, which suggests that 
the electrons excluded from the origin by strong interlayer Coulomb repulsion are 
trapped by the correlation hole in gii(r) within r/i ~ 4. Since the intra-layer 
correlation gn for d/i > 1.6 is almost the same to that of composite-fermion liquid 
state, Agi2{r) represents the correlation of composite fermions between the layers. 
The almost same amplitude of Agi2{r) at r/i ~ 1 and 3 for d/i > 1.6 actually shows 
that the electrons in the other layer bind holes to form composite fermions. 

With decreasing d/i from infinity, the correlations of composite fermions in 
different layers monotonically increases down to d/i ~ 1.2 as shown in the enhance 
of |Z\(7i2(r)| at r/i ~ 1 and 3. But further decrease in d/i broadens the peak at 
r/i 3 in Agi2{r) with the decrease in the correlation hole in Agii{r), and the peak 
at r/i ~ 3 in Agi2{r) finally disappears. This change in the correlation function 
shows how the composite-fermion liquid state evolves into excitonic state: The large 
correlation hole in the same layer, which is a characteristic feature of the composite 
fermions, is transfered into the other layer to form excitonic state. The correlation 
functions in Figs. [16] and [T7] are continuously modified with the decrease in d/i from 
oo to 0, which supports continuous transition from the compressible liquid state to 
the excitonic state. Fig. [TUj also shows that the peak in Agi2{r) at r/i ^ 3 made 
by the binding of an electron to the hole around the origin gradually disappears 
with decreasing d/i from 1.2. This means the gradual break down of the concept of 
composite fermions. 

The break down of the composite fermions around d/i ~ 1.2 affects the character 
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Fig. 18. The change in correlation lunction through the excitation from the ground state to the 
first excited state, i/ — 1 and A'^^ = 18 with Ly/Lx — 1.0. 

of the lowest excitations, which is clearly shown in the level crossing of the excited 
state at d/i ~ 1.2. The change in the character of excitation is confirmed by the 
correlation functions in the excited state. Figure [TH] shows the difference in the pair 
correlation functions gij{r) between the ground state and first excited state defined 



where gfj{r) and g^{r) are the pair correlation functions in the ground state and the 
first excited state, respectively. Sgij{r) in Fig. [18] show that there is a discontinuous 
transition between d/£ = 1.2 and 1.3, which supports the level crossing in the first 
excited state. 

Below d/i ~ 1.2, 6g{r) have large amplitude at r/i ~ 2 and 6, which shows elec- 
trons are transfered between the inside of r/.^ ~ 4 and its outside. Small singularity 
at r/i ^ 5.5 is due to finite size effects of square unit cell. Above d/i ~ 1.2, only 
^gi2if') have large amplitude at r/i ~ 1 and 2, which shows the electrons within 
r/i ^ 4 in different layers are responsible for the lowest excitation. This result sug- 
gests that the low energy excitations are made by composite fermions in different 
layers for d/i > 1.2. 



In this paper we have reviewed the ground state and low energy excitations of the 
quantum Hall systems studied by the DMRG method. We have applied the DMRG 
method to two dimensional quantum systems in magnetic field by using a mapping on 
to an effective one-dimensional lattice model. Since the Coulomb interaction between 
the electrons is long-range, all the electrons in the system interact with each other. 
This fact seems to severely reduce the accuracy of the DMRG calculations. However, 
in the magnetic field, one-particle wave functions are localized within the magnetic 
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length £, and the overlap of the one-particle wave functions exponentially decreases 
with increasing the distance between the two electrons. This means the quantum 
fluctuations are restricted to short-range and the effective Hamiltonian is suited for 
the DMRG scheme. This is the reason why relatively small number of keeping states 
is enough for quantum Hall systems compared with usual two dimensional systems. 

In quantum Hall systems, filling v of Landau levels is determined hy u = Ne/N^, 
where A^,^ is the number of flux quanta and related to the magnetic field as Nff) = 
{e/h)LxLyB. Thus so many types of the ground state are realized only by chang- 
ing the uniform magnetic field B. Since the ground state of free electrons in par- 
tially filled Landau level has macroscopic degeneracy, Coulomb interaction drasti- 
cally changes the wave function. The character of the ground state is sensitive to 
the Landau level index N and the filling u, which modify the effective interaction 
and the mean distance between the electrons. This is the source of many interesting 
low temperature properties of quantum Hall systems and their inherent difficulties. 
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